Percolation Transitions in Scale-Free Networks under Achlioptas Process 
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It has been recently shown that the percolation transition is discontinuous in Erdos-Renyi networks 
and square lattices in two dimensions under the Achlioptas Process (AP). Here, we show that when 
the structure is highly heterogeneous as in scale-free networks, a discontinuous transition does not 
always occur: a continuous transition is also possible depending on the degree distribution of the 
scale-free network. This originates from the competition between the AP that discourages the 
formation of a giant component and the existence of hubs that encourages it. We also estimate the 
value of the characteristic degree exponent that separates the two transition types. 

PACS numbers: 64.60.ah,64.60.aq,68.35.Rh 



The Achlioptas process (AP) is a network evolution 
process in which the number of vertices is fixed as N, 
and edges are added one by one at each time step follow- 
ing a given rule that prevents the formation of a target 
pattern. Recently, Achlioptas et al. [H studied the perco- 
lation transition (PT) for the Erdos-Renyi (ER) model Q 
following an AP rule, called the product rule (PR) in 
which the formation of a giant component is discouraged. 
In their study, the network was developed by choosing be- 
tween one of two randomly selected edges; the selected 
edge had a lower value of the product of the size of the 
two components that edge is joining. They found that 
the giant component emerged suddenly at a percolation 
threshold p c , and that the PT was first-order. This tran- 
sition pattern differs drastically from the continuous PT 
occurring in the conventional ER model. The transition 
is delayed as p c = L c /N w 0.88, larger than p c = 1/2 
for the conventional ER model, where L c is the number 
of edges added to the system up to the transition point. 
More recently, Ziff Q found the same first-order tran- 
sition in the two-dimensional bond percolation clusters 
under AP. Similar explosive transition pattern has also 
been observed in a jamming transition model of Internet 
packets 0). 

Here, we study the PT in a model scale-free (SF) net- 
work under the AP rule. SF networks contain hetero- 
geneous degrees, and their distribution follows a power 
law, Pd{k) ~ k . To construct artificial SF networks,a 
stochastic model called the Chung and Lu (CL) model la 
is used. Similar to the ER model and the static model [6| , 
the CL model starts with a fixed number of TV vertices in- 
dexed i = 1, . . . , N. Then a vertex i is assigned a weight 
of Wi = (i + i — 1)~ M , where (X G [0, 1) is a control param- 
eter, and i oc N 1 ' 1 / 2 ^ Q for 1/2 < fi < 1 and i = 1 for 
fi < 1/2. Then, two different vertices are selected 

with their probabilities equal to the normalized weights, 
Wi I s }2i k Wk and Wj / J2k w k ; respectively, and an edge is 
added between them unless one already exists. This pro- 
cess is repeated until pN edges are created in the system. 



The obtained network is SF in degree distribution with 
the exponent A = 1 + l//i. Henceforth, we will use the 
CL model to study the percolation transition of scale free 
networks in PR (SFPR). 

The mechanism by which a giant component in PT 
forms in conventional SF networks with 2 < A < 3 is dif- 
ferent from that in ER networks. In a ER network, as the 
number of edges L = pN increases in the system, mul- 
tiple isolated small components are created and merged 
together. This process continues up to the finite percola- 
tion threshold p c where a single giant component emerges 
through an abrupt coalescence of those small compo- 
nents. On the contrary, in SF networks with 2 < A < 3, 
the percolation threshold is zero in the thermodynamic 
limit. Thus, the giant component initially develops with 
the largest degree vertex as the seed, and grows continu- 
ously by aggregating small-size components. The devel- 
opment and growth of the giant component is the result 
of relatively high probability of a vertex being chosen 
in the giant component 0. In the SFPR, on the other 
hand, two vertex pairs are selected according to the afore- 
mentioned weights. During the network growth, if two 
vertices get selected from the same component, an edge 
is created between them with no change in component 
size. Thus, the existence of a giant component implies 
that even under AP, the probability of growing the giant 
component is very high. This leads us to ask the fol- 
lowing question: what is the impact of introducing the 
AP rule on the nature of the percolation transition in SF 
networks? 

We obtain the following results by performing exten- 
sive numerical simulations for the SFPR model: There 
exists a tricritical point A c , estimated to be between 
2.3 < A c < 2.4, such that when 2 < A < A c , the tran- 
sition point p c is zero in the thermodynamic limit, and 
the PT is second-order as in conventional SF networks. 
When A > A c , however, p c is finite, and the transition 
is first-order. The jump in the giant component size at 
the first-order transition point decreases as p c decreases. 
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The phase diagram is depicted in Fig. [T] In finite-size 
systems, however, p c {N) is finite even when A < A c and 
the transition is first-order. In addition to this new fea- 
ture, many other unexpected behaviors emerge. 
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FIG. 1: (Color online) Phase diagram of the percolation tran- 
sition in the SFPR network. Here, p — L/N is the edge den- 
sity, and A is the control parameter corresponding to the de- 
gree exponent of non-PR SF networks. A second-order (first- 
order) PT is represented by a solid line (dashed line). The 
tricritical point is denoted as "TP." 



For each case, the edge added to the system is selected 
as follows: In (i), the edge that minimizes the product 
of the component sizes on each side of respective edge 
(PR) is chosen. In (ii), the edge ei is chosen, leading to 
no change in component size. For (iii), an edge is chosen 
randomly between the two. Schematic picture of the se- 
lection rule in the AP is depicted in Fig. [2] Henceforth 
A = l + l//iisa control parameter of simulation: We find 
that in SFPR A is not the resulting degree exponent, un- 
like in the conventional CL model (see below). 
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FIG. 3: (Color online) The fraction G of the giant component 
versus the edge density p for (a) the CL model under AP, 
and (b) the conventional CL model. Data were obtained for 
networks with various control parameters A (2.2, 2.4, 2.6, 2.8, 
3.0, and 4.0 from left (blue) to right (green)). System size 
was fixed at N = 10 7 . 







4> -ft 






(iii-a) 








(i) 




(iii-b) 



FIG. 2: (Color online) A schematic diagram of the selection 
rules in AP for cases (i)-(iii) defined in the text. In case (i), 
two inter-component edges are drawn at random, and one of 
them is chosen to be connected according to the product rule 
(PR). In case (ii), one edge is inter-component and the other 
intra-component edge. The latter is chosen. In cases (iii-a) 
and (iii-b), two intra-component edges are drawn, and one is 
randomly chosen to be connected. 

Specifically, numerical simulations are performed for 
the CL model with the PR. At each time step, two candi- 
date edges, e\ and e-i are drawn from the system with re- 
spective probabilities as described previously and added 
is that that minimizes the product of the component sizes 
on each end of the respective edge. Depending on the 
type of the edge, there are three possible cases: (i) both 
edges ei and ei are inter-component ones, (ii) one edge 
ei is intra-component, and the other ei inter-component, 
or (iii) both edges e\ and ei are intra-component. Two 
subcases of (iii) are shown in (iii-a) and (iii-b) of Fig. [5] 



We measure the fraction of vertices in the giant com- 
ponent, denoted as G, averaged over 10 2 ~ 10 4 different 
network configurations, as a function of p. We define the 
PT point, denoted by p c (N), in a system of finite size 
N as the point at which the local slope of G is maxi- 
mal. This position is consistent with the peak position 
of the susceptibility defined below. We also define the 
discontinuity of G [§], denoted as 5G, as the height of 
the intersection point of two tangent lines, one from the 
rapidly increasing transition region and the other from 
the smoothly increasing curve after the jump. Indeed, 
G shows the first-order phase transition at p c (N) in fi- 
nite size systems as shown in Fig. [3] As the parameter 
A — ► 2 (equivalently p, — ► 1), the transition point p c (N) 
and the jump SG decrease. To understand the behavior 
of G{p) in the N — ► oo limit, numerical simulations are 
performed for various system sizes in Fig. [4j We find that 
there exists a critical value A c , estimated to be between 
2.3 < A < 2.4, such that for A < A c , p c {N) decreases 
to zero as N increases (FigfJJa)) in a power-law manner 
Pc{N) ~ N _1/c with C > (inset of FigH^a)), and thus 
p c {N — > oo) — ► 0. The exponent £ depends on A. For 
example, 0.15(1) for A = 2.2. The jump SG at 

p c (N) decreases to zero as 8G ~ N~@>^, where the ex- 
ponent (3 also depends on A. For example, /3/C ~ 0.23(1) 
for A = 2.2 (FigSpD)). Thus, we conclude that the PT is 
continuous in the thermodynamic limit, and Achlioptas 
suppression is not effective in this case. When A > A c , 
however, p c (N — > oo) converges to a finite value (inset 
of FigfJJd)). The estimated values of p c (N) for different 
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Ns and p c (oo) are listed in Tabic I. In finite size systems, 
p c (N) — p c (oo) ~ N^ 1 ^. For example, estimated value 
of the exponent 1/C = 0.29(1) for A = 2.8. C + 1 in- 
dicates that the first-order transition for A = 2.8 is not 
critical ■ To check the nature of the PT in the thermo- 
dynamic limit, we denote Lq = pqN and L\ = p\N as the 
number of edges at which the value of G reaches 1 / s/N 
and 0.3, respectively. We find that there exist a scaled 
quantity A/N - 8 with A = L\ — Lq, which converges to a 
finite value as N — * oo for A = 2.8 (FigfJJe)). The scaling 
factor N - 8 < N indicates that the transition is of first 
order [l| ■ It is interesting to note that the susceptibility, 
defined as x = XL s2n s with n s , the number of s-size 
components per node and the sum excluding the largest 
component, diverges as TV — > oo even when the transition 
is first-order. We find that x max = x(Pc(N)) ~ 
with 7/C 0.4 and 0.7 for A = 2.2 and 2.8, respectively, 
shown in Figs. |4jc) and (f). Interestingly, j/( w 0.7 
remains unchanged for A = 4.0 and 00. 

Since the second-order and the first-order transitions 
meet at A c , A c is a tricritical point. To estimate the po- 
sition of A c , we measure successive slopes of the function 
p c {N) with respect to N for several values of A and plot 
them as a function of 1/N in the inset of Figj5] We find 
that the successive slopes decrease to zero for A = 2.4 
and 2.5, while they converge to a finite negative value for 
A = 2.3. Thus, we conclude that the tricritical point is 
between 2.3 < A c < 2.4, shown in Fig. [5] 

TABLE I: Estimated percolation threshold p c values for finite 
(iVi = 10 6 and N2 = 10 7 ) and infinite system sizes, and the 
obtained degree exponents A' at p c (oo) for various A. Errors 
in the last decimal points are given in parentheses. 



A 


Pc(JVi) 


Pc(N 2 ) 


Pc(oo) 


A'(p c ) 


2.2 


0.33(1) 


0.23(1) 





2.8(1) 


2.3 


0.42(1) 


0.33(1) 





3.0(1) 


2.4 


0.49(1) 


0.42(1) 


0.18(1) 


3.1(1) 


2.6 


0.60(1) 


0.57(1) 


0.52(1) 


3.5(1) 


2.8 


0.68(1) 


0.66(1) 


0.65(1) 


3.8(1) 


3.0 


0.73(1) 


0.73(1) 


0.72(1) 


4.2(1) 


4.0 


0.83(1) 


0.83(1) 


0.83(7) 


6.3(1) 



The relative frequencies of occurrence of the three cases 
of (i)-(iii) of Fig. [2] during the evolution is related to the 
degree effectiveness of AP. We find that the case (i) oc- 
curs dominantly with a probability nearly one during the 
period p < p c (N), in which an attached edge connects 
two isolated components, merging them into a larger 
component. Above p c {N), it decays rapidly since a gi- 
ant component is already there. The cases (ii) and (iii) 
begin to occur when p is close to p c . Next, we examine 
the component-size distribution during the evolution. In 
early time regime p <C p c (N), the component size distri- 
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FIG. 4: (Color online) (a) Same plot as Fig. 3 but for various 
system sizes N — 10 5 ,10 6 and 10 7 , from right to left. Con- 
trol parameter A = 2.2. Inset: Plot of p c (N) versus 1/N. 
The solid line is a guideline with slope 0.15, indicating that 
p c (oo) — ► 0. (b) Plot of the jump SG around p c (N) versus 
1/N. Solid line is a guideline with slope 0.23. (c) Suscepti- 
bility versus p. Inset: The peak value versus 1/7V. (d) Same 
as (a) for A=2.8. Inset: Same as the inset of (a) for A=2.8. 
Solid line is a guideline with slope 0.0, indicating that p c (oo) 
is finite, (e) Scaling plot of A/N - 8 versus N for A = 2.8, 
where A/N = pi —po with p\ and po being the edge densities 
when the fractions of the giant component reach G = 0.3 and 
G = A r ~ 1/ ' 2 for the first time, respectively, (f) Same as (c) 
for A = 2.8. Error bars in each data point are within symbol 
sizes. 

bution exhibits an exponential decaying behavior. As p is 
increased, the distribution develops a hump in large-size 
region, which is made through the coalescence of small- 
size components, resulting in the abundance of large-size 
components. As p just passes p c , these components fi- 
nally condense into one giant component, resulting in 
the disappearance of the hump and a power-law distri- 
bution of component sizes. This behavior proves that the 
self-organization process operates during the very short 
transition period even under the action of AP. 

Unexpectedly, the tricritical point A c is located in the 
range (2.3,2.4). To understand the underlying mecha- 
nism, we measured the degree distribution of the SFPR 
network. The degree distribution follows a power law 
Pd{k) ~ k~ x . Yet, the exponent A' varies, depending on 
the edge density p for a given A, as shown in Fig. [5J We 
find that A' decreases as p increases. Numerical values of 
A' obtained at p c (N = 10 7 ) as a function of A are listed in 



4 




VN 

FIG. 5: (Color online) (a) Plot of p c (N) versus 1/N for A = 
2.3 (o), 2.4 (o), and 2.5 (□). Error bars in each data point 
are within symbol sizes. Inset: plot of successive slopes of 
p c (N) versus 1/N. For A = 2.4 (o, green) and 2.5 (□, red), 
the successive slopes approach zero, indicating that p c (oo) is 
finite. For A = 2.3 (o, blue), the successive slopes approach a 
finite negative value, indicating p c (oo) = 0. 
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FIG. 6: (Color online) (a) Plot of the degree distribution 
P d (k) of the SFPR network when A = 2.2. We find that the 
degree exponent A' w 2.8 when p = 0.234 ss p c (N) (o, red), 
and A' fs 2.25 when p = 0.8 (o, purple). The degree distri- 
bution of the conventional CL network at p = 0.234 is drawn 
(□, blue) for comparison. The system size is fixed at N = 10 7 . 
Inset: plot of the degree distributions of the SFPR network 
at p c {N) for various system sizes N = 10 5 , 10 6 , and 10 7 from 
left (purple) to right (red), showing the distributions' insen- 
sitivity to system size. Data have been shifted vertically for 
clear view in the main panel, and for easier comparison in the 
inset. 

Table I. Since the degree exponent A' obtained at p c (N) 
turns out to be insensitive to system size N (the inset of 
Fig©, A' at p c (N = 10 7 ) may be regarded as the one at 
p c (oo), even though A' is not defined at p = 0. Interest- 
ingly, when A < A c , A' < 3. Thus, we can assume that 
p c (oo) = when A' < 3 at p c as long as A < A c . This re- 
sult is reminiscent of the well-known fact that p c (oo) = 
when A < 3 in conventional uncorrclated SF networks. 



In summary, we have studied the percolation transi- 
tion in the evolution of SF networks governed by AP. 
The nature of the phase transition changes from contin- 
uous to discontinuous as the degree-exponent parameter 
A is tuned past a tricritical value A c (Fig. Q}. This phe- 
nomenon originates from a competition between AP that 
discourages the formation of a giant component and the 
existence of hubs in SF networks that encourages it. 

Note added in proof: Shortly after the submission of 
this manuscript, we became aware of a similar work [p3T | 
under preparation. It uses a different model from ours, 
the configuration model, exhibiting similar properties 
with some differences. 

This work is supported by KOSEF grant Accelera- 
tion Research (CNRC) (Grant No.R17-2007-073-01001- 
0), and NAP of KRCF. Thank Dr. Fortunato for sending 
us their paper and an anonymous referee for introducing 
Ref. [nj]. 
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